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1.0 INTRODUCTION 

Prediction of .^^design^^comand! communication and 

SS‘ systems° n for “octet vehi f f^^in^fup^ovLde the 
acceptable degree of assurance t certainty of correct command 

required quality of , d ^ a means if the system is to be 

execution must be Two methods have normally ^ een t £° 

certified as reliable. ^ “J™tical analysis, based on the 

perform this functions (1) ^ ntenna element types and their 

known properties of ba j nrt-ive structures of simp 

behavio? in the presence of conductive st models of the 

ScS-s s. -sst 

design process. elec tromagnetic 

The rigorous mathematical applied to only a few 

scattering 1 processes has been successfully appl^ offered 

simple geometric shapes, and h Dro bi e ms that are presented by 

r feasible means for solving the problems influence the 

the large and complex struct ch vehicle or a spacecraft, 

pattern of an antenna mounted on a has relied heavily on 

SESSSSS' performed Sn antenna pattern ra nges.^The ^arr^er 

r h ic a ie te “sl i^chls^ require a very hig^ scaling Uc or, 

fabrication of ^"become very difficult because 

their full-scale counterparts m y 0 ting frequencies have 

■f t.v.pi -r extremely small size. P . -receiving technology 

vehicles f^fTo 

enable ^“^ 0 ^ Sc^an?e^s? 

make 6 the ^r^if fmSd^f Vienna impractical. 

Fortunately , recent years have brought considerable progress 
in two fields related to this problem: 

in aU S ti-i £i ~^ 

SCatt rtmputer technology has progressed at a P-nomenal^ 
The combination of these two pjoduled by Complex 

StluctilS by°mlthlmlticIi e clmputation, and at a reasonable cos 
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The purpose of this task is to identify those computational 
codes which will enable antenna design personnel of the George C. 
Marshall Space Flight Center to perform antenna pattern 
computations for several applications involving the use of 
antennas on physically-complex structures and to facilitate the 
use of these codes by the incorporation of such supplementary 
software as may be necessary. 

Two computer codes have been identified as being best-suited 
to the current needs of the MSFC at this time. These codes have 
been examined and tested and are recommended for use by MSFC. 
Both codes were developed at the Ohio State University. One is 
based on the use of the Geometrical Theory of Diffraction and the 
other is based on use of the Method of Moments. They are 
considered to be complementary, inasmuch as the GTD code is 
suitable for use in problems involving electrically-large 
scattering structures, while the MM code is useful for 
electrically-small scatterers. 


2.0 SURVEY OF EXISTING CODES 

A survey was performed to identify such scattering codes as 
might be suitable for performing the type of computations 
required in the antenna design work at the Marshall Space Flight 
Center. Discussions were held with several personnel who are 
authoritative in the field of electromagnetic scattering theory. 
These discussions led to a general conclusion that most computer 
codes that have been written to solve the problems of interest 
fall into one or more of the following categories: 

1. It is company-proprietary and thus not available. 

2. Little or no documentation is available for other users. 

(It is used by the people who wrote it) 

3. It is written to solve only certain types of problems. 

4. It is not yet ready for release. 

Inquiries regarding specific codes confirmed the above 
conclusion. 

The two recommended codes were found to have the following 
advantages : 

1. The codes were readily available to qualified 
applicants . 

2. They are general enough to be used for a wide variety of 
applications . 

3. They are exceptionally well documented. 
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3.0 


CAPABILITIES AND LIMITATIONS OF THE SELECTED CODES 

The two recommended codes, NEC-BSC and ESP3, are based on 
two entirely different analytical processes. Thus, the ^ r 
strengths and weaknesses are quite different. The NEC-BSC code 
is especially suited to the solution of scattering problems which 
involve electrically- large structures and in which the scattering 
bodies can be approximated as assemblies of fl at plates and 
elliptic cylinders. However, it is not well suited to the 
treatment of electrically-small scattering bodies or scatterers 
which are separated by very small distances as measured in 
wavelengths. Also, it lacks the ability to treat certain 
scattering mechanisms and includes only flat plates and elliptic 
cylinders as basic model shapes. The ESP3 code, on the o er 
hand, is well-suited to solve scattering problems which involve 
electrically-small bodies. Scattering centers which are 
electrically close together do not present any difficulty. Its 
method of solution enbodies all scattering mechanisms. Its cruet 
limitation is its large memory requirement and long run times. 
These features are related to the size and geometric resolution 
of the scattering body (or bodies), so that use of the code is 
practical only for scattering bodies of about three wavelengths 
across or smaller . 


4.0 DESCRIPTIONS OF THE SELECTED CODES 


4.1 NEC-BSC 


This code was developed at the Ohio State University and was 
released for use in December, 1982. It is based on the use o 
the Uniform Theory of Diffraction, an augmented version of the 
Geometrical Theory of Diffraction first introduced by J.B. Kelie 
of New York University in the 1950s. Keller s original techniqu 
offered a simple and efficient means for mathematically 
representing the individual scattering mechanisms that 
collectively comprise the scattering properties of a body, 
enabled the user to model the scattering body as a set of basic 
geometric shapes and to apply the mathematical rela tions betwe 
the incident fields and the scattered fields that had been 
derived for each basic shape. The totalscattered field could 
then be found as the complex sum of the fields scattered by all 
of the basic elements. The GTD method of analysis had two 
significant advantages: 


(1) The amount of mathematical computation required was 
relatively modest; thus a computer code based on this 
method would not impose high demands on the machine to 
be used, nor would it result in long run times. 

(2) The use of this method provides some insight into the 
nature of the scattering process taking place in a 
particular problem. This advantage results from the 
fact that the scattering process is divided into a set 
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of individual scattering mechanisms, each of the 
mechanisms being related to some geometric property of 
the body. The effects of the individual mechanisms 
may be examined separately. 

Despite its attractive advantages, the GTD possessed some 
deficiencies which affected its validity under certain 
conditions. Probably the most important of these deficiencies 
related to the so-called shadow boundaries that were associated 
with diffraction by the edge of a conducting body. The GTD is 
based on a ray-optics concept which treats electromagnetic waves 
as rays, or straight-line segments which emanate from one point 
and which may terminate on another point. When the ray-optics 
model is applied to the conducting edge mentioned above, two 
angles exist at which the value of the diffracted ray abruptly 
changes. It assumes a double value at each of the angles. 
Clearly, this result does not represent reality. 

The shadow-boundary problem was solved several years ago by 
R.G. Kouyoumj ian and P.H. Pathek of Ohio State University. They 
derived a transition function which, when combined with the 
original Keller expressions for calculating diffraction, resulted 
in an equation which produces a diffraction value that varies 
smoothly through the transition regions. It has been shown that 
this function properly represents the diffraction mechanism that 
operates in the transition region. 

Other refinements have been added to the GTD over the years. 
Expressions have been derived to represent diffraction mechanisms 
associated with curved edges and curved surfaces, including the 
contributions of creeping waves. The modern version of the GTD, 
incorporating these improvements , is known as the Uniform Theory 
of Diffraction, or UTD . 

Many, but not all, of the known diffraction mechanisms are 
incorporated in the NEC-BSC code. Those which are included are 
associated with diffraction by the surfaces and edges of finite- 
length elliptic cylinders, flat plates and the wedges formed by 
the joining of flat plates. The scattering mechanisms which are 
computed in the NEC-BSC code are: 

o Singly reflected fields 

o Doubly reflected fields 

o Singly diffracted fields 

o Ref lected-dif f racted fields 

o Dif f racted-ref lected fields 

Double diffraction is not computed in the NEC-BSC code. 
However, warnings are given in the printed output data, giving 
the angles at which double diffraction would occur , and 
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identifying the mechanism which produces it. 

The limited number of basic shapes (two) which the code 
makes available to geometrically describe a scattering structure 
somewhat restricts the ability of a user to faithfully represent 
a vehicle structure. The addition of a cone model and either a 
hemisphere or ellipsoid model would add appreciably to its 
applicability to actual flight vehicles. Perhaps a more serious 
deficiency is the fact that interaction (reflection and 
diffraction) between cylinders and flat plates (the two basic 
shapes) is not included in the computation. Blocking, or 
shadowing of the rays by these bodies is included * k 

Interactions between two cylinders is computed, but only tor a 
special case: that in which the axes of the c ^?; inde T r ?, 

parallel, and only for the plane perpendicular to the cylind 

axes . 

The NEC-BSC code offers a great deal of flexibility to the 
user. Some of the options are: 

1. It allows computation of either near fields or far 
fields . 

2. Backscatter, bistatic scatter or antenna patterns may be 
computed . 

3. Either great-circle cuts or conical cuts may be made. 

4. The pattern coordinate system may be oriented in any 
desired way to the reference system. 

5. Either electric or magnetic source types may be used. 

6. The source current distribution may be selected from bhe 
basic options offered or may be input from the user s tabulated 

data. 

7. Field values may be computed either as a function of 
angle or as a function of frequency. 

The insight which the NEC-BSC code provides regarding the 
scattering mechanisms which operate in a particular 
considered to be a very useful feature. This insight is pr 
by two methods : 

1. The fact that the model of the scattering body is made 
up of discrete geometric elements enables the user to vary ® 
parameters of a particular element, or perhaps to e^inate that 
element, to determine the nature and magnitude of the effect of 
that particular element on the overall pattern. 

2. A selectable feature of the code permits the 
contributions of the individual scattering mec anis 
(reflections and diffractions) to be printed out separately. 
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Thus, the antenna design engineer is assisted in anylyzing the 
scattering process. 

In summary, the NEC-BSC code provides a very useful 
capability to rapidly compute antenna patterns which result from 
the use of various antenna types in complex scattering 
situations, but must be used with care to assure that the output 
data are valid. Caution is recommended in two areas: 

1. Formation of the model should be performed in such a way 
as to achieve the best electrical approximation to the actual 
structure. The model which best describes the total physical 
structure may not be the best model. 

2. The presence or absence of various scattering mechanisms 
in various geometric situations should be considered when 
interpreting the output data. 

4.2 ESP3 

The ESP3 code is based on use of the method of moments. In 
using this method, the scattering surfaces are divided into small 
areas, or "patches". Thus, the computer code is often referred 
to as a "patch code". The scattering structure is assumed to be 
illuminated by a specified incident wave or excited by an applied 
voltage in the antenna structure. An unknown value of current 
flows in each of the patch areas as a result of the incident wave 
or applied voltage. The direction and complex value of this 
current is determined by the direction, amplitude and phase of 
the incident field (scattering problem) or the location and value 
of the impressed voltage (antenna problem) and by compliance with 
the required boundary conditions. When the total field at a 
specified point is set equal to the sum of the contributions of 
the many (unknown) patch currents, an equation in many unknown 
quantities is obtained. If a second point is specified, a second 
equation is obtained. When the number of equations is equal to 
the number of unknowns, the resulting matrix may be inverted and 
the values of the unknown currents may be determined, as well as 
the value of the distant (or near) field which results from these 
currents . 

The ESP3 code may be used for either antenna pattern 
computation or for computing the scattered fields which result 
from a body being illuminated by an incident wave. The 
scattering body (and the antenna structure, for the antenna 
problem) must be expressed by the user in terms of multi-sided 
flat plates and straight wire segments. The flat plates are 
divided by the code into the patch elements described above, and 
the currents which flow in these patches are expressed (as 
dictated by the user) as either surface current density modes or 
as filamentary currents. Currents in the wire segments are 
expressed, of course, as filamentary currents. 

The user defines the structure by specifying the positions 
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of the corners of the plates being used and the positions of the 
ends of the wire segments . Locations of load impedances and 
generators in the structure are also specified by the user. 
Plates and wires may be joined in any variety of ways to 
approximate the actual structure. 

The ESP3 code treats the entire body as a single unit, 
rather than as a collection of separate scatterers, as with the 
NEC-BSC code. This treatment results in an advantage and a 
disadvantage. The advantage is that the computation does not 
depend on consideration of various scattering mechanisms 

associated with separate scattering elements (some of which may 
not be possible to compute). The disadvantage is that the 
insight gained by viewing the contributions of the individual 
scattering elements is lost. The net gain is a more rigorous 
treatment of the scattering problem. 

As compared with the NEC-BSC code, the number of 
computations required when using the ESP3 code is enormous, and 
the computer run time reflects that fact. It is typically forty 
to sixty times as long as that required by NEC-BSC for the same 
problem. A larger working memory is also required. These 
characteristics are largely dependent on the electrical size of 
the model being used, the patch size selected and the choice of 
whether the patch current density mode or the filamentary current 
mode is used. In any case, the model should not exceed more than 
about three wavelengths across in size. 

The ESP3 code should be most useful in two applications: 

(1) Computation of antenna patterns or scattering 
properties of electrically-small antennas and/or scatterers. 

(2) Computation of current distributions on antenna 
elements for inclusion in the input data of the NEC-BSC code. 

Detailed information regarding the ESP3 code may be found in 
the User's Manual for that code. 


5.0 SUPPLEMENTARY INPUT DATA PROGRAMS FOR NEC-BSC AND ESP3 

Two additional codes, INSCAT and INDAT4 have been written by 
Applied Research to provide a more convenient means for entering 
input information required by the two OSU codes. When these 
codes are used, the required input information is requested from 
the user by questions presented on the monitor screen. These 
questions are answered, one at a time, from the keyboard. When 
all questions have been answered, the body of data is stored in a 
data file in the proper sequence and format, ready to be accessed 
by the main program. The use of these codes is optional. 

The data files may also be prepared directly by the user by 
entering the required data in the input data file. The 
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disadvantage of using this method is that the user must take care 
to place all data entries in their proper positions and in the 
proper format in the data file, and must not forget any required 
entry. No cues or directions are provided to the user by either 
OSU program (except by reference to the user's manual) in the 
input process. 

Two command files have also been prepared for each of the 
OSU codes. For the NEC-BSC code, the user may initiate program 
execution by typing the command @ SC ATI . This command calls up 
the SCAT1 command file which performs the following functions: 

(1) It automatically calls the data input program INSCAT 
which requests the input information from the user. 


(2) It assigns input and output files. 

(3) It calls the NEC-BSC executable program SCAT, which 
runs and computes all required output data. 


(4) It deassigns input and output files. 

The user may choose to run the computational code 
without using the interactive input program. this case, tne 

initiating command is 0SCAT2. It is assumed that the input data 
file has already been prepared. The SCAT2 command file performs 
the following functions: 


(1) It assigns input and output files. 

(2) It calls executable code SCAT, which runs and computes 
all required output data. 


(3) It deassigns input and output files. 

For the ESP 3 code, two command files, PATCH1 and PATCH2, are 
also prepared. They are initiated by the commands JPATCHl and 
@PATCH2 , and they are used in the same manner as the 
@ SC ATI and 0SCAT2 are used for the NEC-BSC code. 


6.0 APPLICATION OF THE CODES TO LAUNCH VEHICLES AND SPACECRAFT 


The NEC-BSC and ESP3 codes may be used in at least three 
ways in the analysis of antenna performance on launch vehicles 


and spacecraft: 


(1) The nature and magnitude of the scattering effects 
various structural elements may be evaluated individua Y 
collectively. 


of 

or 


( 2 ) 

antenna 


That portion of a scale model actually needed for 
pattern range measurement may be determined ana y ica y. 
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(3) Complete antenna radiation patterns for antenna/ vehicle 
configurations may be computed, subject to the limitations of the 
code being used. 

The first use mentioned may be highly beneficial when 
decisions regarding antenna type and location are being made. 
Effects on the composite antenna pattern by one or more 
structural elements will be determined in general by the primary 
pattern and polarization of the antenna being considered for use 
and by its location and orientation on the vehicle body. 
Evaluation of these effects by use of a computer should provide 
appreciable savings by reducing the time and cost involved in 
model preparation and pattern range measurement. 

The second use mentioned above will enable the user to 
determine just how much of a vehicle scale model is really needed 
to obtain valid pattern data. In some cases, particularly those 
cases in which a very large vehicle is being scaled and a high 
frequency is being used, a large part of the model may not 
contribute appreciably to the pattern. If the entire body of 
such a vehicle is simulated by a model which has been scaled by a 
factor which enables it to be physically accommodated on the 
range, fabrication of scale-model antennas may be made difficult 
or impossible. However, if it has been determined that only a 
known fraction of the vehicle body contributes substantially to 
the pattern, then only that part of the body may be modeled, so 
that the scale factor does not need to be as great as would 
otherwise be required. Validation should be performed by 
comparison of computed patterns and measured patterns . 

It is stressed that the computer codes produce their 
greatest benefit and least risk when used in conjunction with 
experimental measurements . They should not be thought of as 
stand-alone methods for the solution of antenna design problems. 

6.1 Candidate Problem: Use of Radar Antenna on Shuttle SRB . 

It has been suggested that one of the selected codes be used 
to analyze the effects of the Space Shuttle orbiter , the external 
tank and the solid-rocket boosters (SRB) on the radiation pattern 
produced by a radar transponder antenna located in a specified 
region on the surface of an SRB as shown in Figure 1. This 
problem has been examined and the analytical capability of the 
two codes has been evaluated for the recommended case. 

6.1.1 Feasibility of using the ESP3 code. 

Use of the ESP3 code is not feasible for scattering bodies 
of more than about three wavelengths in size. The shuttle 

configuration has dimensions of hundreds of wavelengths. 

Therefore, this code is clearly not suitable for use in solution 
of this problem. 
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Area considered 
for antenna location 





6.1.2 Feasibility of using the NEC-BSC code. 

The NEC-BSC code, being based on the ^^^rcompStltionl 
diffraction, is best suited to handling scatter ing^comp^ ^ 

which involve structures of the sis ze bei mode ls is the 

has the advantage that one of 1 external tank, both having a 

xs£ s 1 ass ^ 

plates^and^cylinders , but°™4 «9uire * large nUmber ° £ **"“ 

elements for satisfactory fidelity. 

Before proceeding with a ■ discuss ion. 

method, the elements of this P^° ttern da ta will be identified, 
the type and quality of computed pattern data wm u 

(1) The region of the SRB being u“h“ 9 foS 

transponder antenna contains surfac g asses _ the affects of 

of circumferential rings. Tne aonxcy 

these rings has not been determine . 

(2) That portion of the ext *^ a ^ a ^ s W ^ C e X te?nai al lS?f itt 

the antenna location region electrical size of the 

characterized by a corrugated s h a P®- frequency, so that the 

corrugation segments is large at the radar frequency, 

surface cannot be modeled as that of a smooth cylinder, 
it be represented adequately as a surface impe 

,3, The Shapes of the forward and aft : ends of the SRB ^nd 
of the external tank cannot be precisely 
deficiency may not be very important. 

(4) interaction (excluding sh^ bodieS 

SRB bodies and the external tank an 

can be computed only in the roll plane. 

(5) If the orbiter body is model ® d * S \ irs^order° n scatter 
plates and elliptic cylinders only first 

mechanisms will be computed for that body. 


( 6 ) 


Effects of surface coatings have not been evaluated. 


in view of the above cons^tio^ ccTOtl^o^ *£ 
radiation pattern of an SRB anbe ^*'. tion ma y seem a very 
effects of the entire structura ~ omn -i e xitv of the configuration, 
questionable task. Indeed, the p nec-BSC code, dictate 

together with certain limitations of are 

that the task be conducted in a. way £ oblem so that computation 
considered in the formulation of the problem, so 

of invalid data can be minimized. 

It appears that maximum benef if f ects b produced^ by specific 
portions 6 of t^ot^sSre, thus gaining an understanding 
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Of the magnitude and nature ^ 

pattern. An advantage of this annro^ ribUt i° n to the — »P°site 

-Je hi 

B 

Pattern measurements^ or' "thos^sections . perform ^ ^xpe^mlntaf 

6- 1.2.1 Modeling of the Proposed Antenna 

of using the NEC-BSC^ode^n s!l?in!L? SS S SSing the feasibility 
to consider the possibilitv ?f ?? 6 SRB ante nna problem was 

considered for use in this an antenna which is Seing 

in Cr th S€ r ^ av f tY'backed slot antenna^ it?* 6 cand f date antenna is 
ln the C band. Unfortunate! ? ??' Ibs °P era ting frequency is 

capability to model any lo^ce ^ ?° de ?° es not include LL 
surface, and the candidate ant-nn bat 1S i° ca ted on a curved 

by this limitation reaches fa/L/^ type * T he problem 

radar antenna problem <->-f ar beyond the scope of the cdd 

most) launch vlhltle tL 8 ° UrSe ' because a large Sumber (if not 
ouryed surfaces. lus it antennas Ire minted on 

de Y ls J<? to model such antennal in^ lmpor J ant th at some method be 
and which will produce valid data/ W&Y that the code wil1 accept 

it was C ?efi?Ld ng tha? Ur f a SaIf°i nted cavit Y- b acked slot antenna 
mounted on a fi 3f 1 a naif-loop antenna of suitahi^ • ' 

distribution, 3 “ d havlng • 

?hi tern and polarization are concerned 316 "™,"* the slot as far as 
the consideration of the ThlS rea lization led to 

theoretical current segment! ° f Usinq a set of 

a way as to produce ged above the surface in such 

half loop, in ordelto minimi? are a PP r °*i-tely those o 
a rectangular "half-i OOD » the . required number of segments 

only three segments tl configuration was selected, z^Slna 
configuration violate! Fl U u re 2a. However, ev2n thil 

space of at least i /a requirement of any GTD code* th^i- 

°f the system. Two of theL?*/*’/- 8 ^ 3 *' between any two elements 
ground surface. This deficlnov S ln Figure 2a a ^ a touching 
elements to the oosiHi was re — ved by raising these two 
realized that mTSSildo.fS? j" FigUre 2b ‘ « mus? *£ 

!!^n gUratl ° n of this shape s!!h reSent ? n actual antenna 
arrang e m ent would not be achilvabl ar \ . elec trical current 

the basic electrical laws w™ 8 ln . pr actice; it would violate 
theoretical source modei The^m^ 1 ^ quite acceptable as a 
making the end elements short!? /J/ 81 1S fu f ther improved by 
emplitude accordingly. The lnc reasing their current 

modeling process. it permits source mode was used in this 

of current and any orientation < ]’ g ^ ni ®^ t . 0 f any arbitrary value 

^ / h ? Se cur rents d o nit intlZaS' iXV1 ^ al cur rent segments? 
magnitude and phase of thecuwS ^ W1 J h each other - Both 
independently. The capabilit? ?! ^ e ach segment is specified 

paointy to independently specify all of 
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these variables provides a very powerful means for producing a 
theoretical model which closely approximates the performance of 
the actual antenna . 

The three-element configuration described above may be used 
to simulate a single slot. If an identicai configuration is 
used and rotated 90 degrees, the second slot of the* ^ossed ®ied 
antenna is simulated. If the first configuration is so arranged 
that all three of its elements lie in the plane normal to the 
axis of the cylinder over which it is located and the second one 
is arranged so that its elements lie in the plane of the cyl 
axis, we have what will henceforth be referred to 
Configuration 1 and Configuration 2, respectively. 

Two sets of antenna patterns were computed. The first set 
was computed by use of the current- segment mode i in conj unction 

with a 24-inch long by 24-inch diameter c J rc Y^^ Yl J:£ e measured 
patterns were computed in an attempt to duplica . 

patterns provided by the manufacturer of the antenna that is 
being considered. The manufacturer's patterns are shown in 

Figure 3. The patterns seem to indicate that the two slots 
driven with equal power levels and 90 phasing. The second set of 
patterns were computed by use of the same antenna con 
as used in the first set, but with a 12-foot diameter cylind ^ 
simulating the SRB body. Patterns were computed for t^ep 
plane and for the roll plane for each of the two antenna 

configurations . 

Computed patterns for the 24 x 24 inch cylinder will be 
discussed first. Only the dominant polarization component wi 
be constderid in each case, the other component being essentially 
zero for the principal planes. The roll plane cu 
configuration 1 is shown in Figure 4 Th % s ^ on ^ f Cr ^ in ^a?I- 

propagation (diffraction) around the surface of th ® ar !field 
diameter cylinder is evident. The amplitude of 

vicinal is down only about 12 dB. from the peak value at 120 

degrees around the cylinder from the a ^ e ^ na * ^ hi pLANE ter LINEAR 
be compared with the pattern entitled ROLL PLANE, LINEAR 

VERTICALLY" from the antenna manufacturer's patterns S titled 
Fiaure 3, or with the outer envelope of the pattern entitiea 
•'ROLL PLANE, ROTATING LINEAR" The degree of agreement cannot be 
determined with certainty because no dB. scale is shown 
manufacturer f s pattern charts* 

The pitch plane pattern for configuration 1 on the 24 x 24 
inch cylinder il shown in Figure 5. The dominant poUr«.ti« 
component in this case is produced by the current eluent that is 
parallel to the cylinder surface. Therefore its radiated signal 
will fall to almost zero in the directions of the cylinder oe s not 
just as that of a longitudinal slot antenna would do. It does n 
abruptly fall to zero in those directions because of the finit 
length of the cylinder (12 inches in each direction) , as seen 
Figure 5 and in the inner envelope of the pattern in 

PLANE, ROTATING LINEAR" of Figure 3. Note that the pattern 
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OF POOR QUAUT* 


MODEL 8 1 1 5 C 


\ 


vmS I 


rotatFna linear typical RADIATION PATTERNS 
ON A 24“ DIA. CYLINDER 24 LONG 
^ ^ ° N (5.65 GHz) 


PITCH PLANE 

rotating linear 


roll plane 
linear vertically 


pitch plane 
linear vertically 


Figure 3 . Measured patterns of crossed-slot antenna. 



Figure 4 . 


Pattern of antenna model, configuration 1, on 24" x 24 
cylinder, roll plane. 
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Figure 5 is rotated 90 degrees relative to that shown in 
Figure 3. The reason for the radiation in the interval between 
the angles of 60 degrees and 120 degrees is not clear. However, 
the level in this region is 23 dB. or more below the peak level, 
so is perhaps not significant. 

The roll-plane pattern for configuration 2 on the 24 x 24 
inch cylinder is shown in Figure 6. The dominant contribution is 
from the current element that is parallel to the cylinder 
surface, and rolls off fairly rapidly near the 90-degree an ^ 270- 
degree angles, just as that of a circumferential slot a ^enna 
would do. It does not fall abruptly to zero at those angles 
because of the finite size and the curvature of the cylindrica 
surface. This pattern should be compared to the inner envelope 

of “hi pattern entitled "ROLL PLANE, ROTATING LINEAR" in 

Figure 3. Note that there is a single value of 9.6 dB. compute 
at the angle of 180 degrees. All adjacent values are below -30 
dB. This is clearly a computational error. The reason for this 
error is not known. 

The pitch-plane pattern for Configuration 2 on the 24 x 24 
inch cylinder is shown in Figure 7. Using this configuration an 
pattern plane, all three current segments in the configuration 
contribute to the radiated signal to generate the dominant 
polarization component. The situation is analogous to t a 
the circumferential slot antenna, which would ave 
omnidirectional pattern in the half -plane if the ?; ang ^\ 
diameter of the cylinder approached infinity. The pattern 
deviates from that shape because of diffraction at the ends 
the cylinder. Comparing the pattern of Figure 7 with the pattern 
entitled "PITCH PLANE, LINEAR VERTICALLY" in Figure 3 (note the 
90-degree rotation) or with the outer envelope of the pattern 
entitled "PITCH PLANE, ROTATING LINEAR", the patterns are seen to 
be very similar except in the region around 90 degrees in 
Figure 7 (180 degrees in Figure 3). The computed pattern of 

Figure 7 goes to an extremely high value at 90 degrees, whil 
the measured pattern of Figure 3 falls almost to zero ^ . 

region (180 degrees). It may be noted that we are operating in 
the plane that contains the 9.6 dB. anomaly seen in Figure 

6. 1.2. 2 Application of Antenna Model to SRB 

After computation of patterns in the principal planes for 
each of the two model antenna configurations on a 24 x Z 4 men 
cylinder, the two configurations representing the two slo ^ a 
the proposed antenna, and after comparing the computed pa 
with measured patterns, the two model configurations were use 
conjunction with a cylindrical body having the same diame ( 
feet) as the SRB. The length of the body was 20 fe J t : . 

length was chosen, not to represent the actual length of Y P 
of the SRB, but to minimize the effects of diffraction a 
cylinder ends. The flat ends of the model cylinder would no 
represent the actual shapes of the forward and aft ends 
SRB very satisfactorily. Again, patterns were compute i 
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cylinder, roll plane. 


L9 


24" x 24 


Figure 7 


Pattern of antenna model, configuration 2, on 
cylinder, pitch plane. 




configS?atioir eS PitCh) f ° r each of the model antenna 

through n The c °mputed patterns are shown in Figures 8 

Fiaure h 8 ro ^~P la ne pattern for Configuration 1 is shown in 
thp „ ' K c ? m P a rison with the pattern shown in Figure 4 reveals 
Sound hp ba f C Shape ,' ■ exce P t that the radiation does not w?ap 
the Sase of ? yl ; n ^ r to nearly the extent that it does in 

ST?J the cylinder^in 

Figured pit °£ -p l ane Pattern for Configuration 1 is shown in 
p .p ■ * • very similar to the pitch“plane pattern nf 

S fl SL°" 24x2 i lnch blinder, shSwn in Wg"e 5? 

degree to 120-^m erence ls the absence of radiation in the 60- 
7. ^ 120 degree range seen in Ficrure 5 Thic; Hi f f 

edges^a? th^Is^t^^ lon^Sr?* dlffraCti "3 

sh™rii 0 Fig^fr mp exIipfia?^I sgffiss^^ 

appears in s Spt th f fc an apparent computational anomaly 

direction 111 an ? ul f r around the 180-degree 

. Th< f validity of the computed data in this region has 

wculS E. rm i n a d ; Otherwise , the pattern shape is th« Shioh 

as thisone d r ° m a clroumfe rential slot in a cylinder such 

Figured 1 Pi T? h '£ d ? n ? patlern J for Configuration 2 is shown in 
7 fo? the '24 v iT ? co T pared with the pattern shown in Figure 
7 24 x 24 inch cylinder. The left-hand sides of the two 

expect frS^ VerY sl ™ ;Llar ' and are the Patterns which one would 
indi of Se sL?t rOU T fe 5 ent f al slot - The strongly-illuminated 
diffrant^H 6 mal1 c Ylr nder (they are close to the source) produce 

patte r ?f1romTt S o Ts^gL^l^^lhS ^r^e 06 - 3 r j PPle L 

aSenSa and thus 1 "" 96 °\ 1 f nder ' tifends aje far^rom thS 

the Pattern ^ ^ eak } Y illuminated. The right-hand side of 

Figure 7 has bLn°n • clearly V nderstood - The apparent anomaly in 
aZ ina n hP S U f 6d ln Para ^ a P h 5. 1.2.1. The split lobe 
FiauJp? 9 ^ h ^ 9 °- de ^ ree region on the right-hand side of 
degrees and 14? P d£ dS t0 the , pai f of lobes seen in Figure 11 at 38 

the wide seDarat-fnn^'f da f fer ence between the two cases is 
cvlindpr ) ration of the two lobes in the SRB case (long 

24-inch long cylinder Wlt nSrr ° W se P aration in th ® of the 

the apSarent 1 ri?a?L„ P h t r ern (F4gure 12 > was computed to confirm 
PP lation between the lobe positions and the cylinder 
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190 180 170 


Figure 8. 


Pattern of antenna model, configuration 
12-foot diameter cylinder, 20 feet long. 


on 

Roll plane. 
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Figure 9. 


Pattern of antenna model , configuration 1, on 12 foot 
diameter cylinder, 20 feet long. Pitch plane. 


► - 
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190 180 170 


Figure 10. 


Pattern of antenna model, configuration 2, on 12 foot 
diameter cylinder, 20 feet long. Roll plane. 
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Figure 11. 


Pattern of antenna model, configuration 2, on 12 foot 
diameter cylinder, 20 feet long. Pitch plane. 
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length. The cylinder length in this case was increased from 20 
feet to 50 feet. The diameter was unchanged. The two lobes in 
the right-hand portion of the pattern are seen to move closer to 
the axial directions, the same direction that they moved when the 
length was changed from 2 feet to 20 feet. 

It may be concluded from a comparison of the patterns that 
these lobes are related to the diffracted radiation from the 
edges at the ends of the cylinder. The phase variation of the 
excitation along one of these edges is related to the distance 
between the edge and the source. Thus, it may be expected that 
the pattern of the diffracted radiation will be a function of the 
cylinder length. However, the phenomenon is not well understood, 
and requires additional study. 


6. 1.2. 3 Consideration of multiple cylindrical bodies 

Three of the four bodies which make up the Space Shuttle 
configuration are right-circular cylinders terminated in ends of 
various shapes. Although it is recognized that the external tank 
has a surface pattern in the antenna location region for which no 
modeling method is known, it was never less considered instructive 
to examine the composite effects of two cylindrical bodies having 
the same diameters and relative positions as the SRB and the 
external tank. 

6. 1.2. 3.1 Two-cylinder configuration using 24-inch diameters. 

Before proceeding to the shuttle simulation, two patterns 
were computed for the case of the simulated antenna on a 24-inch 
diameter cylinder, with another 24-inch diameter cylinder located 
nearby. The axes of the cylinders are parallel. Each cylinder 

is 200 inches long. The antenna (configuration 1) is located 
midway between the ends of the cylinder above which it is 
mounted. The operating frequency is 5.65 gHz. The computed 
patterns reveal some of the characteristics of the NEC-BSC code. 

The configuration chosen for the first pattern is shown in 
Figure 13. The second cylinder is located directly above the 
antenna, the cylinder axes being 36 inches apart (12 inches 
between the cylinders). The pattern is taken in the roll plane. 

The single and multiple scattering effects produced by the 
two cylinders are very evident in the pattern of Figure 13. The 
value of relative gain rapidly changes at four points around the 
pattern, labeled A, B, C and D. The reason for these changes may 
be seen by inspection of Figure 14. The ray path OCE consists of 
a direct ray OC from the current element to the tangent point on 
cylinder 2, thence around the circumf erence , as a creeping wave, 
to a second tangent point from which it radiates toward E. The 
mechanism consists of a single diffraction, and is computed for 
the entire shadow region shown on the sketch. 
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model, configuration 



The ray path OABD includes a second scattering mechanism. 
The direct ray proceeds from O to a reflection point A on 
cylinder 1, thence to a tangent point on cylinder 2. From this 
point it diffracts around the cylinder to the second tangent 
point from which it is radiated toward D. Both reflection and 
diffraction are involved, so that in this case we have the two- 
step reflection-diffraction mechanism. 

For smaller values of cp , multiple reflections between the 
two cylinders, followed by diffraction from cylinder 1 would 
occur. Since the code can accommodate only a two-step process, a 
reflection-reflection-diffraction mechanism would not be 
computed, nor would any higher order process. The absence of the 
contributions of these mechanisms appears as an abrupt drop in 
the computed pattern amplitude at the edge of the shadow zone. 
This effect can be seen at angles A and B in Figure 13. 

The amplitude changes observed at points C and D in 
Figure 13 occur for the same reasons. Figure 15 depicts the ray 
paths involved at these angles. Path OBDF consists of a 
reflection at B followed by a diffraction at D, producing the 
final ray toward F (reflection-diffraction). However, the path 
OACEG involves two reflections followed by a diffraction 
(reflection-reflection-diffraction), a three-step process which 
the code cannot handle. This contribution is not computed in the 
shadow zone of Figure 15, and its abrupt loss at angle C in the 
pattern (Figure 13) produces a sudden change in amplitude. 

A different arrangement was selected to produce the pattern 
shown in Figure 16. Inspection of that pattern reveals the same 
type of phenomena as those discussed for the previous 
arrangement . 

6. 1.2. 3. 2 Two-cylinder configuration using SRB and ET diameters. 

Two patterns were computed for this two-cylinder case. The 
only difference between the conditions used for computing the two 
patterns is the relative position of the antenna in the plane 
normal to the two cylinder axes. 

Two interesting features are seen in Figure 17 (A). The 
pattern amplitude drops sharply at about 66 degrees. This drop 
is caused by shadowing by the external tank, as shown in 
Figure 18. A second feature appears at 306 degrees. An 
interference pattern suddenly begins at that angle and continues 
through the 306-360 degree region. The mechanism that produces 
the interference may be seen in Figure 18. The following 
scattering effects occur in the vicinity of the 306-degree angle: 

(1) The direct ray from the antenna is seen both above and 
below this angle. 

(2) The direct ray OA is reflected at A for angles up to 
that shown in Figure 18. For larger angles, the direct ray is 
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fleeted ray paths for SRB-ET configuration. 


shadowed by the SRB cylinder body, but the diffracted ray 
produced by its travel around the cylinder is then r 
the ET. This is a two-step process and is computed. 

( 3 ) When ray OA exceeds the angle shown in the 

multiple-scattering combinations come j^e/b^the code The 
more than the two-step process accommodated by the coae. ine 
contributions of these mechanisms are absent m the computation. 

The cessation of the interference process occurs when 
conditions arS such that contributions from the ET no longer 
combine with the direct ray to produce an interference pattern. 

The scattering mechanisms that produce the pattern of 
Figure ?7 (B) are shown in Figure 19. The same Phenomena may be 
seen occuring in this case as were observed for the case of 
Figure 17 (A)? but involving different angles and amplitudes. 


7 . 0 CONCLUSIONS 

A review of existing computer codes for . c ^ ul |^Secraft 

5SES. --SKf I C Sii° h SPa=e 

Flight Center . ^Thele^codes were obtained from the 
tt^ -i Tjar q i Thev are entitled NEC-bbO ana 

code is Y based on the the Uniform Theory of Dif ^h^NEC-BSC 

refinement of the Geometrical Theory Structures 

code is well suited to problems involving scatteri g is 

that are large in terms of wavelengths, whl i? * he ^heSwo codes 
best suited to use with small scatterers. ' 

are complementary . 

Certain deficiencies have been identified in th« ut ion B of 

£SlJS SLtS^S be encountered*"!^ thf ^antenna design 
personnel . 

The potential use of the NEC-BSC code in solving an ha * nt J™* 
design problem associated with the the code can be very 

examined and assessed. It is conclud scattering mechanisms 

n^efullv employed for evaluating the scattering . . 

involved in that problem, and to a limited degree may be used^to 
compute the patterns produced by the compi 
configuration. 

Both of the selected codes should prove highly nt onna 

p^r^rn^ 

fodfr'rni^itudy « 

antenna design study will prove to be useful tools in the desig 
process . 
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Figure 19. Reflected ray paths for SRB-ET configuration. 


APPENDIX A 


DATA INPUT CODES 
& 

COMMAND FILES 
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COMMAND FILES 


SCAT 1 


SCAT 2 


PATCH 1 


PATCH 2 


ORiOtNA'-- tv.'-': ; 

OF POOR QUALITY 


h ASSIGN/ USER SYS*COMMAND SYbSlNPUl 
% RUN INSCAT 

5 ASSIGN F0R005. DAT SYS$INPUT 
* ASSIGN OUTPUT. DAT SYS*OU I PUT 
i RUN SCAT 

5 DEASSIGN SYS$OUTPUT 
5 PEASSIGN SYS* INPUT 


i ASSIGN F0R005. DAT SYS*INPUT 
i ASSIGN OUTPUT. DAT SYS*OUTPUT 
$ RUN SCAT 

% DEASSIGN SYS*OUTPUT 
$ DEASSIGN SYS* INPUT 


* ASSIGN/USER SYS*COMMAND SYS*INPUT 

* RUN INDAT4 

$ ASSIGN FOROIO. DAT SYS*INPUT 
$ ASSIGN OUTDAT. DAT SYS*OUTPUT 

* RUN ESP3 

* DEASSIGN OUTDAT. DAT SYS*OUTPUT 
$ DEASSIGN FOROIO. DAT SYS*INPUT 


i ASSIGN FOROIO. DAT SYS*INPUT 
t ASSIGN OUTDAT. DAT SYS*OUTPUT 
■f. RUN ESP 3 

F DEASSIGN OUTDAT. DAT £YS$GUTPUT 
F DEASSIGN FOROIO. DAT 9YS*INPUT 
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rj o o Ci c* t i n o o o o o o 


PROGRAM INSCAT 


ORIGINAL jPAGti iS 

OF POOR QUALITY 


nATA TK.S2 R i^ PR0VIDES AN INTERACTIVE MEANS FOR ENTERING INPUT 
I HE 0HI ° STATE UNIVERSITY NEC-QSC PROGRAM. THE 

Tn ^n^ STED PARAMETER VALUES ARE WRITTEN TO A DATA FILE F0R005 DAT 
TO BE CALLED BY PROGRAM "SCAT". rUK005. DAT 


THIS PROGRAM WAS WRITTEN BY: 


J. WARREN HARPER 
APPLIED RESEARCH, INC. 
5025 BRADFORD BLVD. 
HUNTSVILLE. ALABAMA 35805 


CHARACTER *3 STRING 
CHARACTER *36 INPUT 
CHARACTER *2 NZERO 

OPEN ( 1, FILE= 'FQR005. DAT ', STATUS= ' NEW ' ) 

STRING = 'CE: ' 

WRITE!*, 1 ) 

1 F0RMAT(1H*. 'Enter Comment (max 36 chars) ') 

READ!*, 2) INPUT 

2 FORMAT (A36) 

WRITE! 1, 3 ) STR I NG, INPUT 

3 FORMAT! A3- A36) 

WRITE!*, 200) 

200 i ^?^ AT( 1X ' ' ENTE R NUMBER INDICATING UNITS TO BE USED -'/15X, '1 
1ETERS ' / 1 5X, '2 = FEET'/15X, '3 = INCHES'/) 

READ!*, 2) INPUT 
STRING='UN: ' 

WR ITE ! 1 , 9 ) STRING 
WRITE! 1,2) INPUT 

WRITE!*, 5) 

5 FORMAT! IX, 'FREQUENCY IN GHZ '• ' ) 

READ!*, 2) INPUT 
STRING='FR: ' 

WRITE! 1,9) STRING 
WRITE! 1,2) INPUT 

7 FORMAT! 12) 

8 FORMAT! II, IX, A1 ) 

9 FORMAT! A3) 


M 


C 

c 

c 


WRITE!*, 10) 

10 FORMAT! IX, ' *** ORIENTATION OF 

STR ING= 'PD. ' 

WRITE!*, 15) 

15 FORMAT! IX, 'THETA, PHI FOR Z AXIS, 
READ!*, 2) INPUT 
WR I TE ! 1 , 9 ) STRING 
WRITE! 1,2) INPUT 


PATTERN AXES ***'/) 


THETA, PHI FOR X AXIS 


( 1 LINE) ' ) 


*** I YVt Uh PAlitHN CUT *** 


210 


WRITE!*, 210) 

FORMAT! IX, 'TYPE OF PATTERN 
1/ ' T = CONICAL CUT 


CUT DESIRED ! GREAT CIRCLE OR CONICAL)'/ 
(CONSTANT THETA)'/' F = GREAT 


non 


ENTER t or f, followed by value of 


2 CIRCLE CUT (CONSTANT PHI)'//' 

3THE FIXED ANGLE, ON ONE LINE. ') 

READ ( *, 2 ) INPUT 
WRITE( 1,2) INPUT 

*** ANGULAR RANGE DESIRED *** 


WRITE(*, 220) „ oattcdw- 

220 FORMAT (IX, ' ANGULAR RANGE DESIRED FOR PATTERN. 

1 ANGLE, ANGULAR INCREMENT') 

READ < *, 2 ) INPUT 
WR ITE (1,2) INPUT 


INITIAL ANGLE, 


C 

C 

c 


c 


C 


C 


*** PLATE GEOMETRY *** 


WR I TE ( *, 40) 

40 FORMAT ( IX 'HOW MANY PLATES?') 
READ ( *, * > NOPLT 
IF (NOPLT. LT. 1) GOTO 82 

STR I NG= ' PG : ' 

DO 80 MP = 1, NOPLT 
WRITE (1,9) STRING 


WRITE ( *, 50 ) MP 
50 FORMAT (IX, 'PLATE NUMBER 
READ ( *» 7 ) NUM 
NZERQ= '0 ' 

WRITE (1,8) NUM, NZERO 


',12,/' HOW MANY CORNERS’?') 


DO 70 ME=1 , NUM 
WRITE ( *, 80 ) ME 

60 FORMAT (IX, 'X,Y, Z POSITION OF CORNER NUMBER 
READ ( *, 2 ) INPUT 
WRITE (1,2) INPUT 
70 CONTINUE 
80 CONTINUE 


12 ) 


FINAL 


*** CYLINDER GEOMETRY *--"* 


C 


82 WRITE!*, 84) 

S4 FORMAT (IX, 'HOW MANY CYLINDERS t ) 
READ(*, *) NCYL 
IF(NCYL. LT. 1 > GOTO 99 
DO 94 N= 1 , NCYL 
STRING='CG: ' 


WR ITE (1,9) STRING 
WRITE!*, 86) N 
86 FORMAT (IX, 'LOCATION 
READ ( *, 2 ) INPUT 
WRITE (1,2) INPUT 


(X, Y, Z) OF ORIGIN, 


CYLINDER NO. 


11 ) 


WRITE (*,88) N 
38 FORMAT (IX, 'CYLINDER NO. 
1ETA-X, PHI-X ' ) 

READ ( *, 2 ) INPUT 
WR I TE (1,2) INPUT 


II,' 0RIENTATI0N'/1X, 'THETA-Z, PHI-Z, TH 


\ - 


on 


WRITE ( *< 90) 

FfiRMAT (IX, 'CYLINDER 


RADII : 


RX, RY'> 


C 


n n n 


READ ( *i 2 > INPUT 
WRITE (1,2) INPUT 

WR ITE ( #. 92) 

92 FORMAT ( 1 Xi 'POSITIONS it ANGLES OF END CAPS'/IX, 'POSITION it ANGLE FO 
1R NEGATIVE END, POSITION it ANGLE FOR POSITIVE END. ') 

READ ( ■*. 2 ) INPUT 
WR ITE (1,2) INPUT 
94 CONTINUE 


*#* SOURCE GEOMETRY *** 


99 WR ITE ( ■*, 100) 

100 FORMAT (IX, ' *** SOURCE GEOMETRY ***'/) 

WR ITE ( *, 110) 

110 FORMAT (IX, 'HOW MANY SOURCE ELEMENTS?') 

READ (•#■,-*■> NW IRES 
STR I NG= ' SG: ' 


C 


C 


C 


C 


DO 180 MS= 1 , NWIRES 
WR ITE (1,9) STRING 
WRITE ( *, 120) MS 

120 FORMAT (IX, 'SOURCE ELEMENT NUMBER', 12/) 

WR I TE ( *, 130) 

130 FORMAT (IX, 'X,Y,Z POSITION OF ELEMENT CENTER?') 

READ ( ■*, 2 ) INPUT 
WRITE (1,2) INPUT 

WRITE ( *, 140) 

140 FORMAT (IX, 'THETA ?-< PHI ANGLES FOR LENGTH i>. WIDTH VECTORS'/' 
1 (THETA-L, PHI-L, THETA-W, PHI-W') 

READ ( *, 2 ) INPUT 
WR I TE (1,2) INPUT 


WRITE ( *» 160) MS 

160 FORMAT (IX, 'TYPE, LENGTH AND WIDTH OF ELEMENT ', 1 2/ ' 

1 -1 UNIFORM CURRENT DISTRIBUTION'/' 

2PIECEWISE SINUSOIDAL DISTRIBUTION'/) 

READ ( *, 2 ) INPUT 
WRITE (1,2) INPUT 

WR ITE ( *, 170) MS 

170 FORMAT( IX, 'EXCITATION (MAGNITUDE, PHASE) hOR tLEMENT ' , 12 ) 
READ ( *, 2 ) INPUT 
WRITE (1,2) INPUT 
180 CONTINUE 

STR ING= ' XQ: ' 

WR I TE (1,9) STRING 
STR ING= ' EN: ' 

WR ITE (1,9) STRING 
CLOSE ( 1 ) 

END 


TYPE: ' / ' 
-2 
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o n o o o no o o o o o o o o o noon 


PROGRAM INDAT4 


THIS PROGRAM PROVIDES AN INTERACTIVE MEANS FOR ENTERING INPUT 
DATA INTO THE OHIO STATE UNIVERSITY ESP3 PROGRAM. THE REQUESTED 
PARAMETER VALUES ARE WRITTEN TO A DATA FILE FOROIO.DAT IN THE 
PROPER ORDER TO BE CALLED BY ESP3. 

THIS PROGRAM WAS WRITTEN BY: 

J. WARREN HARPER 
APPLIED RESEARCH, INC. 

5025 BRADFORD AVE. 

HUNTSVILLE, ALABAMA 35805 

DIMENSION NCNRS(IO) 

DIMENSION X( 10), Y(10), Z( 10) 

DIMENSION NPLA ( 10), BDSK ( 10), VGA ( 10), ZLDA( 10) 

DIMENSION PCN ( 3, 10, 10), IA(10), IB(10), SEGM ( 10) 

DIMENSION IREC < 10) , I PNC 10) , IGSC 10) , NAS C 10) 

DIMENSION NSAC 10) , VGC 10) , ZLDC 10) , IFMC 10) 

DIMENSION IABFPC 10), IABAPC 10) 

DIMENSION VLGC 10) , ZL( 10) 

CHARACTER *1 USEWRS, GEOPRT, DFLT, FARZN, USEPLT, MCOUPL, WRTIMP 
CHARACTER *1 CHFREQ, CHGPAT, CHGRAD, CHGCON 

COMPLEX VLG 
COMPLEX ZL 
COMPLEX VGA 
COMPLEX ZLDA 

WRITE ( 5, 10) 

10 FORMAT ( ' THIS PROGRAM ( ESP3 ) IS BASED ON THE METHOD OF MOMENTS') 
WRITE (5, 20) 

20 FORMAT ( ' 4*#****innnnHHf********* ' ) 

WRITEC 5, 30) 

30 FORMAT ( ' ' ) 

WRITEC 5, 2000) 

2000 FORMAT (IX, 'DO YOU WISH TO CHANGE EXISTING INPUT FILE? Y OR N') 

READ (5, 50 ) CHGFIL 

IF (CHGFIL. EQ. 'N OR. CHGFIL. EQ. 'n ' ) GOTO 8 
ICHG=1 

READ IN VALUES FROM EXISTING DATA FILE 

READ (10, *) NGO, NPRINT, NRUNS, NWGS, IWR, IWRZT, INT, INTP, INTD, INWR, IRGM 
1, IFIL 

READC10, *) IFE, IPFE, NDFE, PHFE 

READ (10, * ) IFA, IPFA, NDFA, THFA 

READC 10, *) ISE, IPSE, NDSE, PHSE, THIN, PHIN 

READ (10, * ) ISA, IPSA, NDSA, THSA 

READ (10, * ) FMC, CMM, A 

IF ( USEPLT. EQ. ' N '. OR . USEPLT. EQ. 'n') GOTO 1530 
READ (10,#) NPLTS 
DO 1520 NPL= 1 , NPLTS 

READ (10, *) NCNRS(NPL), SEGM(NPL), IREC(NPL), IPN(IPL), IGS(NPL) 

DO 1510 NCNR=1, NCNRS(NPL) 

READC 10, *) PCN ( 1, NCNR, NPL ) , PCN ( 2, NCNR , NPL ) , PCN ( 3, NCNR, NPL) 

1510 CONTINUE 
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non 


1520 CONTINUE 

1530 IF ( USEWRS. EG. 'N '. OR. USEWRS. EQ. 'n ' > GOTO 810 
READ (10, * ) IWRZM, IRDZM 
READ ( 10, *> NM. NP, NAT, NFPT , NFS1, NFS2 

C 

DO 1550 1=1, NP 

READ! 10, *> X ( I ) , Y< I ) , Z< I ) 

1550 CONTINUE 

C 

DO 1560 1 = 1, NM 

READ( 10, *) IA< I ) , IB ( I ) 

1560 CONTINUE 


OkmSHiv'. 

OF. POOR QUAUTY 


DO 1570 1=1, NFPT 

READ ( 10, *) I FM ( I ) . IABFP ( I > , VLG ( I ) , ZL ( I ) 
1570 CONTINUE 


C 

DO 1580 1 = 1, NAT 

READ ( 10,*) NAS ( I ) , IABAP < I ) , NPLA ( I ) , VGA ( I ) , ZLDA ( I ) , BDSK ! I ) 
1580 CONTINUE 
REWIND 10 


END OF READING OF DATA FILE 


WRITE!*, 2001 ) 

2001 FORMAT! IX, 'CHANGE FREQUENCY? Y OR N'> 

READ!*, 50) CHFREQ 

IF! CHFREQ. EQ. 'N OR CHFREQ. EQ. 'n ' > GOTO 2002 
CALL RCHNGE ! FMC , 0. , 0. , 1, ICHG) 

READ!*,*) FMC 
WRITE!*, 3200) 

3200 FORMAT!/) 

C 

2002 IF ! USEWRS. EQ. ' N OR . USEWRS. EQ. 'n ' ) GOTO 2006 
WRITE!*, 2003) 

2003 FORMAT! IX, 'CHANGE WIRE RADIUS? Y OR N') 
READ!*, 50) CHGRAD 

IF!CHGRAD. EQ. 'N '. OR. CHGRAD. EQ. 'n') GOTO 2004 
CALL RCHNGE ! A, 0,0., 1, ICHG) 

READ!*,*) A 
WRITE!*, 3200) 


2004 WRITE!*, 2005) 

2005 FORMAT! IX, 'CHANGE WIRE CONDUCTIVITY? Y OR N') 
READ!*, 50) CHGCON 

IF!CHGCON. EQ. 'N '. OR. CHGCON. EQ. 'n') GOTO 2006 
CALL RCHNGE!CMM, 0. , 0. , 1, ICHG) 

READ!*,*) CMM 
WRITE!*, 3200) 


C 

2006 WRITE!*, 2007) 

2007 FORMAT! IX, 'CHANGE PATTERN DATA? Y OR N') 

READ!*, 50) CHGPAT 

IF!CHGPAT. EQ. 'N ' . OR. CHGPAT. EQ. 'n ' ) GOTO 2009 

WRITE!*, 2012) _ n .... 

2012 FORMAT! IX, 'SPECIFY TYPE OF PATTERN CUT'/iCX, ' 1 UKfcA I -c 1 kllE lUi 
1 VAR I ABLE THETA) '/10X, '2 CONICAL CUT (VARIABLE PHI)'/) 

CALL ICHNGE! ICUT, 0, 0, 1, ICHG) 

READ!*,*) ICUT 
WRITE!*, 3200) 
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WRITE ( *» 2013) 

2013 FORMAT (IX, 'VALUE OF CONSTANT ANGLE?') 

CALL RCHNGE ( CONANG, 0. , 0. , 1, ICHG) 

READ ( *, * ) CONANG 

WR ITE ( *, 2014) 

2014 FORMAT (IX, 'ANGLE INCREMENT?') 

CALL RCHNGE( ANGINC, 0. , 0. , 1, ICHG) 

READ ( *, * ) ANGINC 

WRITE!*, 3200) 

C 

2009 WRITE! 5, 2010) 

2010 FORMAT (IX, 'CHANGE PLATE DATA? Y OR N') 

READ (5, 50) MODPLT 

IF ( MODPLT. EQ. ' N '. OR . MODPLT. EQ. 'n'> GOTO 2055 
WRITE (5, 2020) 

2020 FORMAT (IX. 'ADD A PLATE? Y OR N') 

READ (5, 50) ADDPLT 

IF ( ADDPLT. EQ. ' N '. OR. ADDPLT. EG. 'n') GOTO 2030 

NPLTS=NPLTS+ 1 

NPL=NPLTS 

CALL PLATE(NPL, NCNRS, SEGM, IREC, IPN, IGS, PCN, ICHG) 

2030 WRITE (5, 2040) 

2040 FORMAT (IX, 'CHANGE A PLATE? Y OR N'> 

READ (5, 50) CHGPLT 

IF (CHGPLT. EQ. 'N '. OR. CHGPLT. EQ. 'n'> GOTO 2055 
WRITE (5, 2050) 

2050 FORMAT (IX, 'WHICH PLATE?') 

READ (5.*) NPL 

CALL PLATE ( NPL, NCNRS, SEGM, IREC, IPN, IGS, PCN, ICHG) 

C 

2055 WRITE! 5, 2060) 

2060 FORMAT (IX, 'CHANGE WIRE DATA? Y OR N ' ) 

READ(5, 50) CHGWIR 

IF (CHGWIR. EQ. 'N '. OR. CHGWIR. EQ. 'n') GOTO 810 
CALL WIRE ( NM, NP, NAT, NFPT , NSF1, NSF2, ICHG) 

GOTO 810 

********* END OF CHANGE SECTION; BEGIN ORIGINAL INPUT ********* 

8 ICHG=0 

WR ITE ( 5, 12) 

12 FORMAT! IX, 'DO YOU WISH TO SPECIFY ANY OF THE FOLLOWING PARAMETERS? 
1 ( Y or N ) ' / ) 

WRITE! 5, 14) 

14 FORMAT (10X, '* NUMBER OF RUNS TO BE MADE'/lOX, '* PRINT OF MODAL CUR 
IRENTS'/lOX, '* PRINT OF IMPEDANCE MATRIX '/10X, '* NUMBER OF SIMPSONS 
2-RULE INTERVALS FOR: '/15X, ' WIRE-TO-WIRE IMPEDANCES '/ 1 5X , 'SURFACE-P 
3ATCH MONOPOLES '/15X, 'DISK MONOPOLES'//) 

READ ( 5, 50 ) DFLT 

IF(DFLT. EQ. ' Y ' . OR . DFLT. EQ. ' y ' ) GOTO 1005 

APPLY DEFAULT VALUES 

C READ 1: 

NRUNS= 1 
NWGS= 1 
I WR=0 
IWRZT=1 
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INT =4 
INTP=6 
INTD=18 
IRGM=1 
GOTO 35 
C 

1005 WR ITE ( 5. 1010) 

1010 FORMAT (/IX, ' HOW MANY RUNS DO YOU WISH TO MAKE? ' ) 

READ (5,*) NRUNS 
WR I TE ( 5, 1020) 

1020 FORMAT (/IX, 'NUMBER OF WIRE GEOMETRIES FOR EACH RUN?') 

READ (5,*) NWGS 
WRITE ( 5, 1025) 

1025 FORMAT (/IX, 'MODAL CURRENT PRINTOUT: '/6X, '0 NO MODAL CURRENT PRINT 
10UT'/6X, '1 MODAL CURRENTS PLUS WIRE/PLATE GEOMETRY PRINTED') 

READ ( 5, * ) I WR 

C 

WRITE( 5, 1030) 

1030 FORMAT (/IX, 'WRITE IMPEDANCE MATRIX IN PRINTED OUTPUT FILE? Y OR N ' 
1 ) 

READ ( 5, 50 ) WRTIMP 

C 

WR ITE( 5, * ) WRTIMP 


IF (WRTIMP. EQ. 'Y '. OR. WRTIMP. EQ. 'y') GOTO 1032 
I WRZT=0 
GOTO 1036 
1032 IWRZT=1 
1036 CONTINUE 
C 

WR ITE ( 5, *) IWRZT 

C 

WR ITE ( 5, 1040) 

1040 FORMAT (/IX, 'NUMBER OF SIMPSONS-RULE INTEGRATION INTERVALS: '//IX, 'F 
10R WIRE-TO-WIRE IMPEDANCES?') 

READ ( 5, * ) INT 
WR I TE ( 5, 1050) 

1050 FORMAT (/IX, 'FOR SURFACE-PATCH MODULES?') 

READ (5,*) INTP 
WR ITE ( 5, 1060) 

1060 FORMAT (/IX, 'FOR DISK MONOPOLES?') 

READ( 5, *) INTD 
WRITE(5, 1070) 

1070 FORMAT (//IX, 'HOW IS THE WIRE GEOMETRY DEFINED? '/ 10X, '0 BY SUBROUT 
1 INE WGEOM'/lOX, '1 DEFINED VIA THE INPUT FILE') 

READ ( 5, * ) IRGM 

C INPUT NON-DEFAULT INFORMATION 

C 

35 WRITE ( 5, 40) 

40 FORMAT (/IX, 'DOES THE MODEL CONTAIN WIRES? (Y or N>') 

READ (5, 50) USEWRS 
WR ITE ( 5, 45) 

45 FORMAT (/IX, 'DOES THE MODEL CONTAIN PLATES? (Y or N) ' ) 

READ(5, 50) USEPLT 

C 

50 FORMAT (Al) 

IF ( USEWRS. EQ. ' Y ' . OR . USEWRS. EQ. 'y') THEN 
INWR= 1 45 
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C 


END IF 


C 


WRITE ( 5, 60) 

60 FORMAT (/IX, 'SPECIFY REQUIRED 
I EOMETR Y ' / 1 OX , '2 PRINT INPUT 
2. '3 PRINT NOTHING') 


OUTPUT: '/10X, 
PARAMETERS & 


'1 PRINT WIRE & PLATE G 
WIRE/PLATE GEOMETRY ' / 1 OX 


READ ( 5, * ) NPR INT 


C 

C 

c 


WRITE! 5, 70) 

70 REaS?5*5^;f? R s? U WANT T ° PRINT INPUT ° ATA BEF0RE COSTING?') 

IF(I NG0=0 EQ ' /y/ 0R - IFIRST E Q . '«') THEN 

ELSE 
NGQ= 1 
END IF 

WR ITE ( 5/ 80 ) 

80 uSARllBLE X THlTA?vIor P ^ °r PATTERN ^T'/IOX. '1 GREAT-CIRCLE CUT 
READ(5 B *f ICUT A 10X ' 2 C0NICAL CUT * VARIABLE PHI)'/) 

WRITE! 5, 90) 

90 ppah?I ' FAR “ Z °NE ANTENNA PATTERN? (Y or N>'> 

READ! 5, 50) FAR ZN 

WRITE ! 5i 110) 

110 FORMAT! /IX, 'VALUE OF CONSTANT ANGLE ') 

READ! 5,*) CONANG 
WRITE ! 5, 130) 

130 FORMAT! /IX, 'ANGLE INCREMENT'?') 

READ (5,*) ANGINC 
WRITE! 5, 120) 

120 R F EaS?5, ( 50)'pl?0uT WANT ^ ° UTPUT F0R PL0TTING ’ (Y 

IF!FARZN. EQ. 'N ' OR. FARZN. EQ. 'n') GOTO 150 

PARAMETERS SET FOR PLANE-WAVE SCATTER CASE 


ISE=0 
IPSE=0 
NDSE=1 
PHSE=1. 0 
THIN=1 . 0 
PHIN=1. O 
ISA=0 
IPSA=0 
NDSA=1 
THSA= 1 


CONSTANT ANGLE VALUE IS SET FOR EITHER PATTERN CUT 

PHFE=CONANG 

THFA=CONANG 

TYPE OF PATTERN CUT IS SET 

IF! I CUT. EQ. 1 ) THEN 
IFE= 1 
IFA=0 


46 


o o o c.i o o d on on o 


ELSE 
I FE=0 
IFA= 1 
END IF 

C 

c OUTPUT DESIGNATION & ANGLE INCREMENT VALUE 

C 

IF ( PLTQUT. EG. 'Y'. OR. PLTOUT. EQ. 'g') THEN 
IPFE= 1 
IPFA=1 
ELSE 

IPFE=0 
IPFA=0 
END IF 
C 

NDFE=ANGINC 
NDFA=ANGI NC 

C 

GOTO 175 

C 

PARAMETER VALUES SET FOR FAR-ZONE RADIATION 
PATTERN FOR BOTH TYPES OF PATTERN CUT 

150 IFE=0 
IPFE=0 
NDFE=1 
PHFE=1 . 0 
IFA=0 
IPFA=0 
NDFA=1 
THFA=1 . 0 

INPUT MATRIX DOES NOT EXIST FOR THIS CASE 

IWRZT=0 

PLATES ARE MODELED AS PERFECT CONDUCTORS 

CMM=— 1 . 0 

"WIRE RADIUS" SET TO . 001 FOR PLATE 

A=. 001 

INFORMATION INPUT FOR PLANE-WAVE SCATTER 

WR ITE ( 5i 160) 

160 FORMAT ( / 1 X« 'SPECIFY TYPE OF SCATTER COMPUTATION: '/10X, '1 
ITER'/lOX, '2 BISTATIC SC ATTER ' / 1 OX, '3 FORWARD SCATTER ' ) 

C 

READ ( 5< * ) I TYPE 
IF ( ICUT. EQ. 1 ) THEN 
ISE= I TYPE 
ISA=0 
ELSE 
ISE=0 
ISA=ITYPE 
END IF 
C 

IF (PLTOUT. EQ. ' Y '. OR . PLTOUT. EQ. 'g'> THEN 
IPSE=1 
IPSA=1 
ELSE 


BACKSCAT 
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IPSE=0 
I PSA=0 
END IF 


C 

C 

c 


r 

C 

C 


C 

C 

C 


— INPUT ANGLE INCREMENTS FuR PLANE-WAVE SCATTER 

NDSE=ANGINC 

NDSA=ANGINC 

INPUT CONSTANT ANGLE FOR PLANE-WAVE SCATTER 

PHSE=CONANG 

THSA=CONANG 


INPUT ANGLE OF INCIDENCE 


WR I TE ( 5. 170) 

170 FORMAT (/IX, 'WHAT IS THE ANGLE OF INCIDENCE? 
READ ( 5, * ) THIN/ PHIN 


(THETA, PHI ) ' ) 


C 

C 

C 


INPUT FREQUENCY 


175 WRITEC 5, ISO) 

ISO FORMAT (IX/ 'FREQUENCY IN MEGAHERTZ?') 
READ (5,*) FMC 


c INPUT WIRE CONDUCTIVITY & RADIUS IF WIRES USED 

C 

IF ( USEWRS. EQ. 'N '. OR. USEWRS. EQ. 'n'> GOTO 210 
C 

WRITE ( 5. 190) 

190 FORMAT (IX, 'WIRE CONDUCTIVITY IN MEGAMHOS/METER?') 

READ (5,*) CMM 
WRITE ( 5, 200) 

200 FORMAT (IX, 'WIRE RADIUS IN METERS?') 

READ (5,*) A 

C 

GOTO 215 

Q 

c SUPPLY FICTITIOUS NUMBERS FOR WIRE PARAMETERS 

C 

210 CMM=- 1 . 0 
A=0. 001 


c 

c 

c 

c 


c 


215 IF (USEPLT. EQ. 'N '. OR. USEPLT. EQ. 'n ' ) GOTO 300 
INPUT PLATE INFORMATION 


WR ITE ( 5, 220) 

220 FORMAT (/IX, 'HOW MANY PLATES?') 

READ ( 5, * ) NPLTS 

WR ITE ( 5, 225) 

225 FORMAT (/IX, 'INPUT TYPE OF PLATE TEST MODE: '/10X. 
1ATCH MODE'/lOX, '1 FILAMENTARY TEST MODE') 

READ ( 5, •»■ ) IFIL 


C 


WR I TE ( 5, 230) 

230 F0RMAT(/1X, 'PLATE INFORMATION — 
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/ ) 


'0 


FULL-SURFACE P 


’ ’ C 



c 

c 

c 


DO 420 NPL=1 , NPLTS 

CALL PLATE ( NPL, NCNRS, SEGM, IREC, IPN, IQS- PCN. ICHG) 

420 CONTINUE 

REUSE OF IMPEDANCE MATRIX 

Disk storage of impedance is not anticipated; therefore: 
300 I WRZM=0 
IRDZM=0 


IF (USEWRS. EQ. ' N ' . OR . USEWRS. EQ. 'n ' ) GOTO 800 
C 

C INPUT WIRE INFORMATION 

CALL WIRE(NM, NP, NAT, NFPT , NSF1, NSF2, ICHG) 

CALL WRPTS( I, NP, X, Y, Z, ICHG) 

CALL ENDPT ( NM, IA, IB, ICHG) 

CALL FPPLT ( NFPT , IFM, IABFP, VLG, ZL, ICHG) 

C 

DO 700 1=1, NAT 

CALL FPNPL ( I, NAS, IABAP, NPLA, VGA, ZLDA, BDSK, ICHG) 

700 CONTINUE 
800 CONTINUE 

C 

c PARAMETER VALUES WRITTEN TO DATA FILE 

C 

810 WRITE( 10, 1 00 ) NGO, NPRINT, NRUNS, NWGS, IWR, IWRZT, I NT, INTP, I NTD, INWR, IR 
1GM, IFIL 

WRITE! 10, *) IFE, IPFE, NDFE, PHFE 
WRITE! 10, *) IFA, IPFA, NDFA, THFA 
WRITE! 10, *) ISE, IPSE, NDSE, PHSE, THIN, PHIN 
WRITE! 10,*) ISA, IPSA, NDSA, THSA 
WRITE! 10, *> FMC, CMM, A 
C 

IF ! USEPLT. EQ. 'N'. OR. USEPLT. EQ. 'n / ) GOTO 760 

C 

WRITE! 10, *) NPLTS 

C 

DO 750 NPL=1, NPLTS 

WRITE! 10, *> NCNRS(NPL), SEGM(NPL), IREC(NPL), IPN!NPL), IGS1NPL) 

DO 725 NCNR=1, NCNRS(NPL) 

WRITE! 10, *) PCN! 1, NCNR , NPL ) , PCN <2, NCNR, NPL ) , PCN! 3, NCNR, NPL ) 

725 CONTINUE 
750 CONTINUE 
GOTO 770 
C 

C INPUT FILLER NUMBERS FOR NO-PLATES CASE 

C 

760 WRITE! 10, *) 0 

WRITE! 10, *) 4, 0. 2, 1, 3, 0 

DO 780 1=1,4 

WRITE! 10, *> 1.0, 1.0, 1.0 

780 CONTINUE 


770 WRITE! 10, *) IWRZM, IRDZM 

C 

C DO NOT WRITE WIRE INFORMATION TO FILE FOR NO-WIRES CASE 

IF ! USEWRS. EQ. ' N '. OR . USEWRS. EG. 'n*') GOTO 1000 
C 

WRITE! 10, *) NM, NP, NAT, NFPT, NFS1, NFS2 


DO 775 1=1, NP 
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WRITE( 10. *) X(I).V<I)»Z(I) 

775 CONTINUE 

C 

DO 111 1=1. NM 

WRITE ( 10. *) I A( I > , IB ( I ) 

111 CONTINUE 

C 

IF ( NFPT. GT. 0 ) GOTO 3000 
WRITEUO,*) 1, 0. ( 1. . 0. ) . ( 1. . 0. ) 

GOTO 3100 

C 

3000 DO 850 1=1. NFPT 

WRITE ( 10, *) IFM< I ) , I ABFP ( I ) , VLG< I ) , ZL ( I ) 

850 CONTINUE 

C 

3100 DO 900 1=1. NAT 

WRITE ( 10. * ) NAS ( I ) , I AB AP < I ) , NPLA ( I ) . VGA( I ) . ZLDA ( I ) 
900 CONTINUE 

C 

100 FORMAT ( 1215) 

1000 CONTINUE 
C 

END 


BDSK ( I ) 
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SUBROUTINE FPNPL(NAT, NAS- IABAP, NPLA- VGA- ZLBA. BDSK, ICHG) 
DIMENSION NAS ( 10), IABAP (10), NPLA ( 10) - VGA ( 10) , ZLDA ( 10), BDSK ( 10) 
IF ( NAT. EQ. 0) GOTO 800 
DO 700 1=1. NAT 
WRITE ( *, 655) 

655 FORMAT (/IX, 'NUMBER OF WIRE SEGMENT ATTACHED TO PLATE:") 

CALL ICHNGE(NAS< I ), 0, 0, 1, ICHG) 

READ ( *, * ) NAS ( I ) 

WR ITE ( *> 750) 

WRITE < #, 660) „ A . tl _ u 

660 FORMAT (/IX. 'ENDPOINT ATTACHED TO PLATE: '/10X, '0 POINT A'/lOX, 

10INT B'> 

CALL ICHNGE( IABAP( I ), 0, 0, 1, ICHG) 

READ (-*»•*) IABAP ( I ) 

WRITER. 750) 

WRITE ( *» 665) 

665 FORMAT (/IX, 'PLATE NUMBER?') 

CALL ICHNGE(NPLA< I ), 0, 0, 1, ICHG) 

READ ( *, * ) NPLA ( I ) 

WRITE ( *, 750) 

WR ITE < *, 670) 

670 FORMAT (/IX, 'COMPLEX 
READ ( *, * ) VGA ( I ) 

WRITE ( *, 675 ) 

675 FORMAT (/IX, 'COMPLEX LOAD IMPEDANCE AT ATTACHMENT POINT / ) 
READ<*, *) ZLDA ( I ) 

URITE<*, 680) 

680 FORMAT (/IX, 'DISK RADIUS 


GENERATOR VOLTAGE AT ATTACHMENT POINT?') 


IN METERS?') 


READ («■,*) BDSK ( I ) 
700 CONTINUE 
750 FORMAT ( / ) 

800 CONTINUE 
RETURN 
END 


c 
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c 

SUBROUTINE FPPLT ( NFPT , IFM, I ABFP; VLG, ZL, ICHG) 

DIMENSION IFM ( 10), I ABFP (10), VLG< 10), ZL( 10) 

COMPLEX VLG 
COMPLEX ZL 

C 

WRITE (5, 560) 

560 FORMAT (/IX. ' INPUT LOCATIONS OF FEED POINTS: '/) 

DO 650 1=1, NFPT 
WRITE (5, 570) I 

570 FORMAT (IX, 'FEED POINT NO. ',13) 

WRITE <5, 575) 

575 FORMAT (IX, 'IN WHICH WIRE SEGMENT?') 

CALL ICHNGE( IFM( I ), 0, 0, 1, ICHG) 

READ ( 5, * ) IFM(I) 

WRITE < ■*, 700) 

WRITE ( 5, 580) 

580 FORMAT (IX, 'AT WHICH END OF SEGMENT? '/ 1 OX . '0 ENDPOINT A'/lOX, '1 E 
1NDP0INT B'/) 

CALL ICHNGE( IABFP( I ), 0, 0, 1, ICHG) 

READ (5,*) I ABFP ( I ) 

WRITE <*, 700) 

C 

WRITE (5, 590) 

590 FORMAT ( / 1 X, 'INPUT COMPLEX VOLTAGE OF GENERATOR AT FEED POINT: ') 
READ < 5,*) VLG ( I ) 

WRITE ( *, 700) 

WRITE ( 5, 600) 

600 FORMAT (/IX, 'INPUT COMPLEX LOAD IMPEDANCE (OHMS): ') 

READ (5,*) ZL ( I ) 

WRITE ( *, 700) 

700 FORMAT ( / ) 

650 CONTINUE 
RETURN 
END 
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SUBROUTINE 
IF( ICHG. EQ. 
IF ( N. EQ. 1 ) 
IF ( N. EQ. 2) 
IF ( N. EQ. 3) 
10 FORMAT (IX, 
20 FORMAT (IX* 
30 FORMAT (IX. 

WRITE ( *, 40 
40 FORMAT (IX, 
50 CONTINUE 
RETURN 
END 


RCHNGE(R1 < R2 
0) GOTO 50 
WRITE(*< 10) 
WRITE ( *> 20) 
WRITE ( #» 30) 

'THE PRESENT 
'THE PRESENT 
'THE PRESENT 
) 

'NEW VALUE (S)?' ) 


, R3, N, ICHG ) 
R 1 

Rl, R2 
Rl, R2, R3 
VALUE IS 
VALUES ARE: 
VALUES ARE: 


FB. 2) 

' /2F8. 2) 
' /3F8. 2) 
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SUBROUTINE ICHNGE(I1, 12, 13, N, ICHG) 

IF ( ICHG. EQ. 0) GOTO 50 
IF(N. EQ. 1 ) WRITE ( *, 10) II 
IF(N. EQ. 2) WRITE**, 20) 11,12 
IF(N. EQ. 3) WRITE**, 30) 11,12,13 
10 FORMAT (IX, 'THE PRESENT VALUE IS ',15) 

20 FORMAT* IX, 'THE PRESENT VALUES ARE: '/2I5> 
30 FORMAT* IX, 'THE PRESENT VALUES ARE- '/3I5) 
WRITE**, 40) 

40 FORMAT (IX, 'NEW VALUE ( S ) ? ' ) 

50 CONTINUE 
RETURN 
END 
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SUBROUTINE PLATE (NPL, NCNRS/ SEGM. IREC, IPN, IGS, PCN, ICHG) 

DIMENSION NCNRS ( 10) , SEGM (10) , IREC (10). IGS(lu). IPN( 1 0 ) . PCN ( 3, 10. 10) 
WRITE (5. 240) NPL 

240 FORMAT (IX. 'PLATE NUMBER 13) 

WRITE ( 5, 250) 

250 FORMAT (IX, 'HOW MANY CORNERS?') 

CALL ICHNGE(NCNRS(NPL) , 0, 0, 1, ICHG) 

READ ( 5, * ) NCNRS ( NPL ) 

WRITE ( *, 272) 

WRITE ( 5, 252) 

252 FORMAT (IX, 'SIZE OF PATCH?') 

CALL RCHNGE ( SEGM ( NPL ) , 0. , 0. , 1 , ICHG ) 


READ ( 5, * ) SEGM ( NPL ) 

WRITE ( *, 272) 

WRITE < 5, 254) 

254 FORMAT (IX, 'IDENTIFY SHAPE OF PLATE 
1ECT ANGULAR ' / ) 

CALL ICHNGE( IREC (NPL) , 0, 0, 1, ICHG) 
READ ( 5, * ) IREC ( NPL ) 

WRITE ( *, 272) 


-'/10X, '0 POLYGONAL '/10X, '1 


*** POLARIZATION SELECTION IS SET TO 3 FOR ALL CASES *** 
IPN ( NPL ) =3 


WRITE (5, 256) 

256 FORMAT (IX, 'WHICH IS THE GENERATING SIDE?') 

CALL ICHNGE ( IGS(NPL) » 0, 0, 1, ICHG) 

READ < 5, * ) IGS(NPL) 

WRITE**, 272) 

WRITE(5, 260) 

260 FORMAT (/IX, 'INPUT POSITION OF EACH CORNER (X, Y, Z> ) 


DO 275 NCNR=1, NCNRS(NPL) 

WRITE ( 5, 270 ) NCNR 

CALL RCHNGE ( PCN ( 1 , NCNR, NPL > » PCN ( 2, NCNR, NPL ) , PCN ( 3, NCNR, NPL ) » 3, ICHG 

READ ( 5, * )PCN( 1 , NCNR, NPL ) , PCN( 2, NCNR. NPL ) , PCN ( 3, NCNR, NPL ) 

WRITE**, 272) 

272 FORMAT*/) 

275 CONTINUE 
RETURN 
END 
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c. 

SUBROUTINE WRPTS( I. NP, X, Y. Z, ICHG) 

DIMENSION X< 10). Y< 10). Z( 10) 

WRITE (5, 520) 

520 FORMAT ( /IX. 'INPUT LOCATIONS (X.Y. Z) OF WIRE POINTS: ') 
DO 530 1=1. NP 
WRITE (5. 525) I 
525 FORMAT (IX, 'POINT ',13) 

CALL RCHNGE ( X ( I > , Y ( I ) . Z ( I ) . 3, ICHG) 

READ < 5,*) X(I), Y ( I ) » Z ( I > 

WRITE < 600) 

530 CONTINUE 
600 FORMAT < / ) 

RETURN 

END 
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430 


440 


450 


460 


SUBROUTINE WIRE(NM, NP, NAT. NFPT, NSF1. NSF2. ICHG) 

WRITE (5. 430) 

FORMAT ( /IX. 'NUMBER OF WIRE SEGMENTS:' ) 

CALL ICHNGE(NM, 0. 0, 1. ICHG) 

READ ( 5. * ) NM 
WRITE (*, 600) 

WRITE ( 5. 440) 

FORMAT (IX. 'TOTAL NUMBER OF 
CALL ICHNGE(NP. O. 0. 1. ICHG) 

READ (5.*) NP 
WRITE ( *, 600) 

WRITE(5, 450) 

F0RMAT(1X. 'NUMBER OF WIRE-TO-PLATE 
CALL ICHNGE(NAT. 0, 0, 1. ICHG) 

READ (5,*) NAT 
WRITE < *i 600) 

FORMAT? lxf ^NUMBER OF FEED POINTS IN THE WIRE STRUCTURE?') 
CALL ICHNGE(NFPT. 0, 0, 1. ICHG) 

READ (5,*) NFPT 
WR ITE ( ** 600) 

WR ITE ( 5. 470) 


POINTS IN THE WIRE STRUCTURE- 


ATTACHMENT POINTS?') 


C 470 FORMAT (IX, 'IS MUTUAL COUPLING TO BE COMPUTED?') 
READ (5, 50) MCOUPL 
50 FORMAT (Al) 

IFtMCOUPL. EQ. ' Y '. OR. MCOUPL. EQ. 'y ) GOTO 480 

NSF1=0 

NSF2=0 

GOTO 510 


C 


480 WRITE (5, 490) 

490 FORMAT (/IX, 'LOCATION OF FIRST FEED PORT. 
CALL ICHNGE(NSF1, 0, 0, 1, ICHG) 

READ ( 5, * ) NFS1 
WR I TE ( •*, 600) 


) 


WRITE (5, 500) 

500 FORMAT (IX, 'LOCATION OF SECOND FEED PORT :' 
CALL ICHNGE(NSF2, 0, 0, 1, ICHG) 

READ (5,*) NFS2 
WR ITE ( *, 600) 

600 FORMAT ( / ) 

510 CONTINUE 
RETURN 
END 


) 
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c 

SUBROUTINE ENDPT<NM. I A, IB, ICHG) 

DIMENSION IA( 10) , IB ( 10) 

535 FORMAT^*/ IX, ' INPUT ENDPOINT A S< ENDPOINT B OF EACH WIRE SEGMENT: '/) 
DO 550 J=1 , NM 
WRITE (5, 540) J 
540 FORMAT (IX, 'SEGMENT ',13) 

CALL ICHNGE ( IA( J), IB( J), 0, 2, ICHG) 

READ< 5, * ) I A ( J ) , IB ( J ) 

WRITE ( *» 600) 

550 CONTINUE 
600 FORMAT ( / ) 

RETURN 

END 
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